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I discuss the relation between harmonic polynomials and invariant 
theory and show that homogeneous, harmonic polynomials correspond 
to ternary forms that are apolar to a base conic (the absolute). The 
calculation of Schlesinger that replaces such a form by a polarised 
binary form is reviewed. 

It is suggested that Sylvester's theorem on the uniqueness of Maxwell's 
pole expression for harmonics is renamed the Clebsch-Sylvester theo- 
rem. 

The relation between certain constructs in invariant theory and an- 
gular momentum theory is enlarged upon and I resurrect the Joos- 
Weinberg matrices. 

Hubert's projection operators are considered and their generalisations 
by Story and Elliott are related to similar, more recent constructions 
in group theory and quantum mechanics, the ternary case being equiv- 
alent to SU(3). 
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1. Introduction. 

There has recently been increased activity in some aspects of classical spherical 
harmonic theory partly in response to the very accurate measurements of the cos- 
mic microwave background and the need to analyse these in a significant fashion. 
One approach has lead to the rediscovery of Maxwell's way of picturing spherical 
harmonics through their 'poles', {e.g. [1-5]), the uniqueness of which description is 
the content of Sylvester's theorem. 

Because of the relation between this topic, and angular momentum theory in 
general, and invariant theory, I thought it might be helpful to set down a personal 
resume of some of these things in a more historical context, if only to draw attention 
to sometimes forgotten, but relevant, work. I therefore try to draw together mostly 
existing material. I also include some peripheral constructions that I feel are worth 
exposing again, such as the Joos- Weinberg matrices. 

The heyday of classical, constructive invariant theory was the second half of 
the 19th century and the first quarter of the 20th. However, it never really went 
away and has recently undergone a resurgence of interest (e.g. Olver, [6]) finding 
applications in coding, computer aided design and pattern recognition.^ 

It is possible to find parallels (and equivalences) between many constructions 
and processes in classical invariant theory and those in 'modern' applied mathe- 
matics. Although my pointing out some of these is probably more of interest than 
utility, I hope to provide something of value. 

2. Spherical harmonics. 

Discussions of spherical harmonics are many. The term itself seems to date 
to the treatise of 1867 by Thomson and Tait, [8], where fairly general definitions 
can be found, but the subject had its beginning in a work of Legendre of 1782 
followed by Laplace's paper of 1782 on potential theory (as we now call it). Later 
fundamental papers are those of Green (1828) and Gauss (1841). Maxwell, [9], lists 
some standard technical references of his time. 

As might be expected, there are a number of approaches. Here I just try to 
give some relevant facts without bothering too much about logical ordering. 

For this reason, I begin with the addition theorem for (solid) spherical har- 
monics which can be regarded as 'classic', being derived, effectively, by Legendre. 

^ Olver's introduction in [6] and Sturmfels' in Hilbert, [7], supply useful historical comments. 
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It is, 

(r'r)^P^(cos7) = C^(r')Cf(r) (1) 

where 7 is the angle between r and r', rr' cos 7 = r . r' . 

I give here, because it's convenient, the useful composition law for these stan- 
dard solid harmonics, Cff{r), in my conventions, 

and the 'inverse', 

^1 ^ ^ ^2 ^ ^ \^Mi M^L^J ^ ' V y ^3 ^ ^ ^ ^ 

The left-hand side of (1) is a homogeneous bipolynomial, of degree 2L, in r 
and r' and is homogeneous of degree L in r and of degree L in r'. It is symmetrical 
in r and r' and is a solid spherical harmonic in either set being sometimes referred 
to as the 6ma;a/ harmonic of r and r', ( [8] ^, [10] ). P^(cos7) is a Laplace coefficient 
(c/ MacMillan, [11] p. 377 ^) and the explicit form of the left-hand side of (1) can 
be found from the classic series for as, (c/ [11] p. 383 equn.(2)), 

(r'r)^P^ (cos 7) = Hj^iv.v')) 

Conversely one could derive this expression from first principles, using rota- 
tional invariance and the harmonic condition, and then deduce the standard series 
expansion of P^, c/ [12]. 

The vectors r and r' can be taken complex, subject to complex orthogonal 
transformations (complex rotations or, equivalently, Lorentz transformations). An 
important case is when one of them, say r', is isotropic, or null, r' = 0, which can 
be achieved by treating the spin-one quantity r', = a, as composed from a 2-spinor, 



iV2 iV2 



^ Thomson and Tait, [8] p. 159, derive (1) by means of Taylor's theorem. See also Hobson, [10] 
§90. 

^As noted by MacMillan on p. 300, (1) can be thought of as a transformation from a harmonic 
with its pole on the ^-axis to one whose pole is on the line {x', y' ,z'). 
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sometimes called the Cartan map, [12]. 

{a^,ay,ag) are the Cartesian components of a and I have chosen a normalisation 
in keeping with the standard/contrastandard definitions. This accounts for the 
differences with [12]. I should therefore say that my convention is that upstairs 
indices correspond to contrastandard, and downstairs ones to standard behaviour, 
in the terminology of Fano and Racah [13], to whom I generally adhere for signs, 
factors of i etc. However my 3j-symbols are the more usual ones of Wigner and 
my angle and rotation conventions follow those of Brink and Satchler, [14]. The 
relation between their surface harmonics and my solid ones is 

1 also give my raising and lowering convention, which is the transpose of Wigner's 
and is the same as that adopted by Williams, [15]. Specifically, for a spin-j object, 
(f), I set 

^ . (6) 

Using (5), (1) reduces to, 

^(a.r)^==Ci,(a)Cf(r). (7) 

C^{a) has the familiar, monomial expression (e.g. Weyl, [16], van der Waerden, 
[17], Wigner, [18], Bargmann, [19], Schwinger, [20]), 

(2L)! ^L+M^L-M 



(8) 



2^/'^! [(L + M)!(L-M)!]^/' 



^ ^(2L)! (^) 

where is (for j half-integral) what I have termed, [21], a null (2j + l)-spinor. 



J — 

The factors here^ are chosen so that the scalar product of with another 
null spinor. 



1/2 



^ Landau and Lifschitz, a very useful book, make the same choice, [22] eqn.(97.4). See also Fano 
and Racah, [13], App.F. 
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is the simple power of a bracket, (c/ [12] (6.150)), 
In terms of this notation, (7) reads. 



a (3 



2j 

(10) 



^(a.r)- = e£)cf(r), (11) 

quoted, with conventional phases, by Schwinger [20]; see [12] 6.149. The left-hand 
side can be considered to be a generating function for the spherical harmonics (c/ 
Erdelyi et al, [23], §11.5.1). Likewise, and rather trivially, the right-hand side of 
(10) can be taken as a generating function for the null spinor ^(-j). Just set t = a/ j3 
and compare with [23] equn 11.7(10). 

If r is also a null vector, r = b, constructed from the 2-spinor (^) by the 
Cartan map, then, from (10) and (11), or directly from (5), 

(a.b)^ = (r7a-e/3)'^ (12) 

As an interlude, I now rederive (11) in order to introduce some related material 
that I feel deserves revisiting. I start from the equation for the product of 2j null 
3-vectors, (see [21] equn. (9)), 

where the t*i - *2j are the spatial components of the Joos- Weinberg matrices, t^^^ = 
f^J■l■■■^^2j which are higher-spin analogues of the space-time Pauli matrices, a'^. They 
are symmetric and traceless on (//). pj is the normalisation, 

p. = 2-^^e-'''^ {4jy}/\ (14) 

They form a complete set (another one!) of (2j + 1) x (2j + 1) hermitian 
matrices and can be expressed as symmetrized powers of the angular momentum 
matrices, examples being, (beware: I use 'j ' in two roles), 

^0...0^;^. ^^O...O^iJ^. fjO...O^^_fhji^jji xijY (15) 

J 2j - 1 Vj J 

The t^-^\ when sandwiched between (2j + l)-spinors, transform like 4-tensors. 
Equation (13) is a higher-spin version of the Cartan map (5). The simplest case is 
j = 1/2, identical to (5), 
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where (Tj are the standard spatial Pauh matrices. 
I now contract (13) with rs to give, 

{a.vf^ =p^l{r,r^...f^-)^. (17) 

From the tensor operator (adjoint) transformation of the t*-' -- one can easily 
show that, [24], 



r r ■ 



Further, the important multiplication law for null spinors is, 

1712 Jl + J2/ 

essentially by Wigner-Eckart. This means, in particular, that, 

(ii\>d'jh'"''')=0, (20) 



Jl 32 "^3 

unless J3 = Jl +J2- 

So, from (19), (17) reads, 

(a. r)^^- = {Aj + 1)V2 e-i {-rf^ Tr {uf^{j) V^^\^v)) , (21) 
where {J) are the matrices discussed by Racah and have matrix elements, 

[<{J)]:' = {Jm \u^{i)\Jm') = [^l^^ . (22) 

The trace on the right-hand side of (21) is a hyperspherical harmonic,^ which 
is a spin-j SU(2) representation matrix expressed in polar coordinates, (x, ^, ??), on 
the 3-sphere. It is evaluated on the equator of S^, x = 7r/2, and is just a surface 
harmonic, C. Generally, 

TV («f V\g)) = i'^Hj^M Cf (e, V) , (23) 
where the radial function is related to Gegenbauer polynomials, 

(2J-L)! ^^/^ 



HjAx) = LI 



(2J + L + 1)! 



(2zsinx)^C2^/i,,(cosx). (24) 



^ These are standard objects and will form the basis of a further work. I refer here now only to 
Talman, [25]. 
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In (21), X = 7r/2, J = j and L = 2j so one needs, 

H,M^/^) = (^^f'mj^-^^- (25) 

Therefore, from (23), using (25) and (14), 

(a . rf' = j^?^ e-^- i2if^ {-rf^ C^.(r) 4^^) 

- o2.7 Pj c'^fy] 

= e— ^■C2-(r)^^2,)_ 

I have thus regained (11), as promised, with integer M = m. 

I wish to prove something that might not be immediately apparent i.e. that 
the quantity on the right-hand side of (13) is traceless, as it ought to be since the a 
are null. The spatial t* '-^ themselves are not traceless but the null nature of ^ saves 
the day as I now show. Firstly, the space-time traceless property gives the spatial 
trace, 

^ii^...i2j_2 _ ^00ii...i2^_2 ^ 

which contains at most a product of {2j — 2) angular momentum matrices and the 
claim, therefore, is that the quantity, 

f'^r^...j'r.^^i) , n<2j-2, (26) 

vanishes. Writing the J* as 3j symbols, the products can be recombined in turn 
using 6j symbols until the ^s are coupled by the same 3j symbol and the null 
theorem (20) can be applied. Because is never greater than 2j — 2, the desired 
result follows. Alternatively, (19) can be applied directly to (26) to give, 

^.TV(u2,.(J)J^i...J^"), 

in terms of the Racah u matrix, (22), which can be shown to vanish. The J are 
spin-j matrices. I leave the details as an exercise. 

I give some further details about the Joos- Weinberg matrices (Weinberg [26]) 
in particular on the relation (18). In fact this relation, or rather the one from 
which I shall shortly derive it, was used by Weinberg to compute the matrices as 
products of angular momentum matrices. A direct algebraic method, which is easily 
implemented is the following. Setting, 

e-^"-J = f;c,(^) (n.J)\ 

fc=0 
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for ^ = TT, only even powers occur for integral j and odd powers for j half-odd 
integral and one has to invert the set of equations, obtained by diagonalisation, 



(-1)^ = J2 ™^''c2fc(7r), m = 1, . . . j , j integral 

fc=0 
J-l/2 

V2 = J2 m^''+'c,,+,{7:), m = 1/2, . . .j , j half - odd integral, 

fe=0 

with CQ(7r) = 1. Torruella, [27], provides an analytic inversion as a function of j but 
it is easier (and equivalent), for any particular j, to use the tabled formulae, 

(_,,m-i/2 ,,n^2n+ M^'- (1/2)') jm'- (3/2)^) . . . ^ - (n - 1/2)^) 

^ ^ ^ (2n+l)! 

71 = ^ ' 

For example, 

j = l, C2(7r) = -2 
3 

j = -, Ci(7r) = 7^/3, C3(7r) = -4i/3 

j = 2, C2(7r) = -8/3, C4(7r) = 2/3. 

From the computed expansion, using (18), the spatial Joos- Weinberg matrices 
can be extracted as sums of symmetrised products of the J\ if desired. I will not 
do so here but refer to the calculation in Weinberg's paper, of which the above is a 
particular case. (Just set = that reference.) 

The t^^^ are actually geared to the relativistic situation. They transform under 
Lorentz transformations, regarded conveniently as complex 3-rotations, as 

g(a-i/3) . J e("+»/3) ■ J = A/ . . . (27) 

where A is the 4-vector Lorentz transformation. The first thing to note is that t^^---^ 
is a scalar under the spatial rotation subgroup, o; = 0, and so can be set equal to 
1 by Schur's lemma. With this in mind, the next step, [26], is to saturate the fi 
indices with the four-vector, q^, and choose A to be a pure Lorentz boost, /3 = 0, 
that kills the 3-vector part, i.e. ^ q'^ — {q, 0), where q' — qo — |qp. Then (27) 
becomes, 

q^ ... tM- = g2j^-2aq.3 ^ ^28) 
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where 

q q 

q = — j- , sinh a = — . 

|q| Q 

Expansion of the right-hand side of (28), as a polynomial in q^, allows the t^^-* to 
be read off, as mentioned before, [26]. 

To obtain (18) from (28), formally continue Qq to zero. Then sinho; = i so 
a = m 12 and we are back to ordinary rotations (through tt!). 

Before leaving this topic, I briefly mention two uses of the t^^\ Firstly, the 
massless equations (Weyl, Maxwell, linearised vacuum gravity, Dirac's higher spin) 
can all be written in the form, [28] , 

t'^-^a,0(a;) = O, 

and, secondly, the quantity ^^.-^ t^^^ 4>{j) is a higher-spin analogue of the Bel- 
Robinson tensor, [24] . 

3. Invariant Theory. 

If ^ and r\ are viewed as variables, and a and /3 as coefficients, the bracket power 
in (10) is the symbolic form of a binary 2j-ic, traditionally written {aiXi ^a^Xr^^ , 
n = 2j {e.g. Grace and Young [29]). For example, for the general quadratic (spin 
one), axf + 2b x^ x^ + cx\, one has a = af, b = a^a^ ^'^d c = a^, but ac 7^ 1? and 
the coefficients, (a, a/2 6, c), do not constitute the (spherical) components of a null 
3-vector (3-spinor), 

Comparing with (5), it cannot escape notice that the Cartan map had al- 
ready appeared algebraically in the 1850's, when Aronhold introduced the symbolic 
method into invariant theory. Since this theory plays a part in what I wish to say, 
I present some expository material. I will not avail myself of the various render- 
ings of classical invariant theory into modern algebraic terms, preferring the earlier 
language. 

The general 2j-ic, ^^j, formed with the coefficients 0^^, is written. 

In the case that the coefficients form a null (2j + l)-spinor, say = «(_,), 
$27 takes the simple power form (10) and the 2j roots of $27 = 0; fo^^ ^ = ^/^; are 
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all equal to a/ (37 The necessary and sufficient condition for all roots to be equal, 
is that the Hessian of the quantic vanish identically. The Hessian for the quantic 
(29) is given in my notation, by 

and so one has the theorem, [21], that if, 

then (^Q) is null. The converse is contained in (20). 

It is fundamental that subjecting the two spinor, '0 = ( ^ ) , to a transformation 
belonging to SL(2,C), and requiring the binary n-ic, (29), to be invariant, 
induces a (particular) Lorentz transformation on the coefficients, (f)^. Invariant 
theory treatments can be found in TurnbuU, [30], Chap. VIII §8, Glenn, [31] §1.2.4, 
[29], §16.^ I repeat some standard details here out of interest. 

The basic transformation i/j ^ ijj' , expressed as the inverse for convenience, is. 

The imposed invariance of $ implies, setting ^' /rj' — t, 



= r7'"$J^',M + mi,M + m2) (31) 
^^'^ e*^-^-$„(^,mi,m2) 

where 1 . is the polarisation operator (or translation operator or directed deriva- 
tive) , 



^ This is an extreme example of a nullform, as defined by Hilbert to be a form all of whose 
invariants vanish. For this one needs only a root of multiplicity [n/2] + 1, n being the form 
order. 

^ The group in invariant theory does not have to be unimodular. For those requiring the basics of 
invariant theory in short, but expert, compass I recommend Dickson's book, [32], which treats 
both symbolic and non-symbolic methods and employs group theory concepts. 
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I interject a small, necessary calculational point. The normalisation of the 
coefficients, 0"^, differs from that of the coefficients as usually defined in invariant 
theory. The two conventions are exhibited in the forms. 



1/2 



m=—j ^ ' 

n 



(33) 



r=0 

with n = 2j and r = j — m, so that, 

</^-^= (^)'^'a,, (34) 

connects the two sets of coefficients. 

The normalisation for the 0s corresponds to the use of the 1-j symbol as the 
raising and lowering metric in weight space, (6), as in (29), and elsewhere, and this 
is an appropriate place to make some elementary remarks on duality which I will 
need later. 

The spinor ip = (—rj,^) is dual to V' = (^)- I can make this projective by 
setting — ?7 —>■ u and ^ — > f , and then refer to tt, v as line, or tangential, coordinates. 
^ and T] are 'point coordinates'. 

The point binary form (33) can be written in line coordinates as follows, 

^2,(0,e,^)= Q!-^^) ' (-1)^— 0-e+-(-r/)^-- 



27 



j-mj (35) 

1/2 



= (-l)2^$2,((/>,«,^;) 



.J+m j-m 



corresponding to (29). The sign factor could be absorbed into the coefficients of the 
line form. 

Putting (33) into the invariance (31) yields 



n 

E(^)a;^'^-^ = e*'-^-*n(a,mi,m2) 



r=0 
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and expansion gives the relation between the old and new coefficients, 



a 



; = -(l.Vj"-'^$Ja,mi,m2), (36) 



I.e. 

ilm 



""^^ (l.Vj^-+-$,,(^,m„m,). (37) 



(2j)!(i + m)! 

If one rather uses rf = t^' , then the alternative formula 
, (n-r)! 



nl 

i.e. 



(m.Vi)'$,(a,Zi,Z2), (38) 



arises ([30], [31]). 

Infinitesimal behaviour and projection operators. 

As mentioned, the induced transformation on the coefficients is that of the 
spin-j representation of the rotation group (when SL(2,C) is reduced, inessentially, 
to SU(2)) and these relations are more or less identical to a standard method of 
finding the representation matrices, V\ of the rotation group, e.g. [33-35,19]. It is 
worth expanding on this in the following way. 

A binary covariant, K, is an SL(2,C), and hence SU(2), invariant, i.e. {cf 
(31)), 

Useful information follows from the infinitesimal expression of this invariance. 

The binary (spin-half) realisation of the generating angular momentum oper- 
ators is, (see Sharp, [36], Bargmann, [19]), 

acting on functions, f{^,r]). A commutation relation is 

[J+,J_] = 2J,. (41) 
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For ease of comparison, I have reverted to the notation J_|_ = ± iJy. These 
generators serve also for the binary invariant group, SL(2). 
The total angular momentum (Casimir operator) is 

and the simultaneous eigenfunctions, /4, of and are easily determined by a 
standard procedure. From (42), they must be homogeneous of degree 2j, where 
j{j + 1) is the eigenvalue of J^, with j a half-integer. The eigenvalue of equals 
m, with —j<m< j. As usual, the can be determined, using the lowering J_, 
from the 'highest weight' function /j, which itself satisfies J+ fj =0 and so is fixed 
to be proportional to the power, . Then, lowering by induction, yields, see (9), 

fj — dii) 
J m Sm ' 

and we see that the are function space variants of the quantum vectors, | jm). 

The matrix elements of in this basis are the usual ones. I do not give them 
because they are contained in the expressions for the transformations induced by 
(40) on the coefficients of the quantic, (29). These already occur in early invariant 
theory in the guise of the Caylcy-Sylvester-Aronhold operators, Q, and O. 

Acting on a binary covariant, of the ground form, (29), (33), the operators 
J_i_ are equivalent to the annihilators, and 0, i.e. , 



(43) 



with, in classic notation. 



n— 1 rj n—1 rj 

fc=0 fc=0 

This theorem was proved by invoking the invariance of the covariant under the 
infinitesimal shears, ^ ^ ^ + Ary, rj rj and ^ ^ ^, J] + ?], generated by J_ 

and J_|_, respectively. (See e.g., Salmon, [37] §148. Dickson, [32,38], gives a succinct 
treatment and Elliott, [39] Chaps. VI, VII a more lengthy one.) The operator 
generates the scalings that maintain ^rj. Olver, in his nice, modern reworking of 
classical invariant theory, [6], discusses the Lie algebra aspects of Q, and O. 
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Transcribed into SU(2) notation via (34), equation (44) reads, {n = 2j), 



1/2 ,m 9 



d(hm-l 
m=l-j 

C= E (0'-™ + l)(^'+™))'^'<^™~'i (45) 

m=l-j 
m=-j 

in which the angular momentum matrix elements can be recognised using the rela- 
tion 

between the Lie derivative operators, X^, and the matrix generators, G^, in the 
(matrix) representation in the carrier space of which is a vector. Thus 



^3 ^ '^0 = e ^3 = [eg"G„ j j _ 

The general structure of the SU(2) covariant, of degree g and order vj — 2s 
is a Clebsch-Gordan coupling of g spin-j objects, to a spin-s quantity (the set of 
coefficients) which is then coupled to a spin-s null spinor (the variables) to a zero 
spin resultant. The invariance conditions (43) are then equivalent to the standard 
invariance of the 3-j symbols under SU(2). 

The operators O and O perform important functions in the construction of 
concomitants and thence in the proof of the Hilbert finite basis theorem. The 
algebra of Q. and O, developed for this purpose, can be translated into angular 
momentum terms using (45) and (46). 

Corresponding to J^, (41), there is the commutator ('alternant') 

-[0.!il = -i:P*-nKa^ = - (47) 

fc=0 m=—j ^ 

originally constructed by Cayley, [40]. Therefore, in total, fl', O' and |[(9',ri'] 
correspond to J + ^cr. 

I reckon weight following Cayley's first scheme, [40], in which has weight 
r — n/2, and ^ and 77 have weights ±1/2. This is more in keeping with the group 
representation value than the one usual in classical texts, being the eigenvalue of 
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Jj,, and has the advantage that the commutator |[0',f2'] is the weight, or scahng, 
operator since (j/^ has weight —m. 

A basic result, essentially due to Hilbert, is that the quantity, 

(a; ,) = (l - ^ + - ^ + . . . ) F(a; ,) , (48) 

where F is a rational, integral, isobaric ^ function of the a^. and ^,r] of zero total 
weight, is a covariant of the ground form (33), or zero.^*^ 
The operator, 7Y_|_, in (48), and its partner, 

take the forms of general projection operators constructed from raising and lowering 
operators so that one could refer to K as the covariant projection of F. It is a 
theorem that every covariant can be obtained in this way. 

The operator equivalent of this projection, customarily couched in angular 
momentum terms, is that the terminating polynomial, 

oo (_-\\r 

^r!(l + r)! + ^ ^ 



is equivalent, using (41), to 



Y -TtH^T J- (51) 
^r!(l + r)! + ^ ' 



when acting on an eigenstate of with zero eigenvalue, m = 0, and projects onto 
the identity, or trivial, representation, i.e. to a scalar. This last statement accords 
with the scalar quality of the covariant, K. The m = condition is just the zero 
weight one on F. 

The operator, (50), is a special case of a projector derived by Lowdin, [41], 
equn.(32) and rederived by Shapiro, [42]. This more general operator is equivalent, 
for quantics, to ( (a)^ is the Pochhammer symbol), 

l!(2m + 2)i ^ 2!(2m + 3)2 3!(2m + 4)3 ^ ' ' 7 ' ^ ' 



^ Isobaric means that each separate term of an expression has the same weight. Olver defines it 

to mean invariance under the scahng subgroup, which is more restrictive. 
10 See Elhott, [39], §182. 
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where, for simplicity, I consider actions on functions of just the coefficients, e.g. 
on homogeneous, isobaric functions i.e. gradients., G{a.), of excess 2m, a number 
which determines the scahng of G via 

I [0,n]G = mG. 

For these terms, consult Elliott, [39]. 

The sufficiency of dealing with just the projection onto the identity has been 
exploited by Shapiro, [43], for arbitrary groups. The basic notion is one commonly 
used in selection rule calculations that computing how much of a given irrep there 
is in, say, a product of two irreps, is equivalent to coupling the three irreps to the 
trivial one. This corresponds to a simple group average, e.g. [35]. It is the same as 
the choice one has in invariant theory of dealing either with the form as a whole (a 
scalar) or with just its coefficients (a carrier space vector). 

As an example of a higher group, Noz and Shapiro, [44] , give the identity pro- 
jection for SU(3) as a polynomial in the generators. In the invariant theory setting 
this involves an extension to the ternary domain, geometrically more interesting, 
but more complicated. The most systematic and elegant treatment of projection 
operators is given by Story, [45] , and a summary in [46] . As an example, the group 
GL(3) is generated by 6 shears (determinant one), and 3 expansions. In Elliott, [39], 
Chap. XVI, the shears are the annihilators, denoted with i,j = x,y, z; i j, 
and the expansions are the non-zero commutators, H[, with ^H'- = 0. 

Elliott effectively works out the structure constants of sl(3) in the Cartan basis, 
the H'- being the Cartan subalgebra in precisely the modern notation, e.g. Racah, 
[47], equn.(72). For example, in terms of the Cell-Mann matrices A3 ~ H'^ and 
Ag ~ H'l — H2 in the fundamental, quark representation, {x,y,z). H[,H2 and H'.^ 
generate scalings of the pairs (y, -2), {z,x) and {x,y) respectively, corresponding to 
the three SU(2) 'binary' subgroups which give F-spin, [/-spin and /-spin. These 
scalings leave the products, yz, zx and xy invariant, each pair as in the binary case, 
leading onto the definition of weights used by Story. 

The projection operators are assembled out of the which are the raising 
and lowering operators, usually denoted by in Lie algebra theory. A comparison 
of the projectors derived by Story, [46], equns.(8), (9) and Elliott, with those of Noz 

Elaboration of this will form part of a further communication. 

For those who wish to quickly refresh their knowledge on the Lie algebra of GL(n)D SL(n), I 
suggest Racah, [47], and Gourdin, [48], esp. Chap. 3. 
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and Shapiro, [44], equns.(2.4), (2.17), shows agreement. See also Asherova and 
Smirnov, [49]. All methods involve chains of subgroups e.g. SU(3)dSU(2). 

There is a close connection between the construction of Lie group/algebra rep- 
resentations and that of covariants from semi-covariants,"*^^ in particular from one, 
given, term of the covariant (its 'source') by repeated action of O' or (9', as in 
Roberts' theorem in the binary domain. (I recommend [32], [39] and [7]). The 
various matrix elements will emerge as in (45) for sl(2). The projection operator 
method for simple Lie groups is given in more detail by Aserova et al, [50]. 

4. Polars and angular momentum addition. 

The polar operators l.Vj^, (32), and their obvious extensions to higher do- 
mains, play important roles in invariant theory. 

The introduction of u and v is not necessary in the binary domain and so is not 
usually made in the classical works. One can easily work directly with (— ?7, Line 
coordinates are needed however for ternary and higher forms and their introduction 
in the binary case allows a more unified treatment, cf Todd, [51], §§1.6, 2.12. For 
example I now bring in the famous 'Omega process' used in the composition of 
invariants, and covariants, to give others. To be a little more systematic, I denote 
point variables by (i = 1, . . . ,p) and the dual line coordinates by m*. A p-ary 
n-ic is written symbolically as, 

/ : = (a.x)" = (o/xj" = {a^x^ + a^x^ + ... + a^x^ 
F : tt^ = (u.a)" = (w^Qij)" = {u^ai + u'^a2 + . . . + tt^ctp)" 

in point and line coordinates respectively. I am most concerned with p — 2 and 
p — 3. 

The geometrical interpretation of forms becomes more involved the higher p is. 
One interpretation-*^^ follows on taking the as homogeneous point coordinates in 
a (p — l)-dimensional space and the vanishing of a form yields a co-dimension 1 
variety (a hypersurface) in this space. In the binary case this is a range of points, 
in the ternary case, a curve and quaternary, a surface. This interpretation provides 
a useful language and visualisation, at least for low p. 

^'^ This will be considered elsewhere. 

^'^ These are the analogues of the extreme weight states, | jj) and | j, —j). 
Not the only possibility. 
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Consider two distinct p-forms of generally different orders, f = = u'^ = F 
and g = b'^ = = G, where I have indicated both the point and line forms. 
Following Clifford, [52], define the generalised polar oi G in respect of /, 
{n-m)\ , _ ™ ^„ _ (n-m)! 



mini ' ^1 ^1 

{n-m)\ , (53) 



= (a./3rar- = a^ar-, n>m. 



The first line is in Grace and Young (for the ternary case), [29], §241, the second 
is Clifford's definition (also in the ternary case). The last line is the symbolic 
expression. In the ternary case, h gives the polar locus, of order (n — m), of G in 
respect of /. 

h is zero if n < m. However in this case, we can reverse the roles of / and G 
and compute the generalised polar of / in respect of G, 

= (!i^(a.v.)%^ (54) 

= (a.;^)^^;;^-" = a;^w^-", m>n, 

which represents the polar envelope, of class (m — n), of / in respect of G.^^ 

For binary forms, the situation is summarised, explicitly in point coordinates, 
by Sturm, [59], who also describes the geometry in terms of ranges of points on 
the complex projective line in a standard way. For this purpose, it is expressive 
to rewrite h, say, using the factored form of tt^, which is, reverting to the original 
binary notation, 

m 
k=l 

SO that setting Vj, — > and Uj, — > —rjj^ and also, according to the definition (53), 
tt — > 5^, V — > d^, 

(n — mV. TTT f ^ d d \ „ 



■"^^ A class curve is often called an envelope or a tangential curve and an order curve is a 'locus' 
A 'curve' can usually be expressed in either of these dual ways. The best known example is 
the conic, which is of order 2, or equivalently, class 2. Simple descriptions, in English, can be 
found in school textbooks such as Robson, [53], Milne, [54], Sommerville, [55]. Possibly the most 
attractive is Askwith, [56]. More advanced is Salmon, [37]. The rather rare book by Scott, [57], 
should be mentioned and the more modern Todd, [51], and Semple and Kneebone, [58]. 
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which shows h as a mixed polar of / with respect to the points (roots) 

Writing h as the form (cx)"^""^, the coefficients, c^., can be determined in 

terms of and b^. They are given in Sturm, for example, [59], and would allow 

one to calculate the 3j symbols ("f "^/^ {n-m)/2yi7 

In the binary case, the quantity is the operator introduced by Cayley and 

which, in point coordinates, is expressed as a determinant^^. Use of line coordinates 

allows it to be written as a contraction, Q = V„ . V^, (sec Todd, [51], §2.12.3). Two 

binary forms, / and g, can be combined into a third using to give their rth 

transvectant, 

(/,ff)-- ^^''^,'^^~'^' ^-(/g)- (56) 
nl ml 

In particular the polars are the 'end values', 

h=if,9r, h' = if,gr, 

and generally {f-,gY is a form of order (m + n — 2r) on making the replacement 
{ui, U2) (^3^2, Xi) (or a form of class {m + n — 2r) on making the dual replace- 
ment {xi,X2) {u2,—Ui). The order (class) runs from (m + n) to \m — n\ and 
transvection corresponds to addition of angular momenta via Clebsch-Gordan, or 
3j, symbols. As a well known example, the second transvectant of a form with itself 
is proportional to its Hessian, H, (30). 

I can now introduce the important notion of apolarity, originally developed in 
the binary case but later extended to higher domains. 

The straight analytical definition is that two forms, / and G, are said ^° to be 
apolar if their polar, h, defined in (53), vanishes. I take n>m for ease. 

An important fact is that the form / is apolar to any form that contains G 
as a factor and whose class does not exceed n. The proofs, which are simple, 
are given in Grace and Young, [29] pp.225, 304. In the binary case this result 
corresponds to elementary angular momentum addition rules, as I now explain. / 
and g correspond, respectively, to spins = n/2 and = m/2 as explained earlier. 

As might be expected, these 3j symbols have simple explicit forms, e.g. Hamermesh, [35], p. 375 
eqn. (9-123). They occur later. 

In the ternary case the operator O, as defined here, occurs on p. 296 of [29] as Q. 

The notion appears to be due to Reye. I will not dilate on its general significance but just refer 
to Coolidge, [60] pp.368, 410, Grace and Young, [29], Semple and Kneebone, [58], Todd, [51]. 

Whether we say apolar to G or to its dual, g, is, metrically, immaterial and I will be somewhat 
slack in making this distinction. 
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The polar h corresponds to spin — ^3 = (n — m)/2. If it vanishes, so that / and 
G are apolar, one has the 3j equivalent 

fm.gm, f Jl k Jl - i2 Q ^ .g^N 

\mi 1712 * / 

i.e. 2{ji — + l = n — m+1 conditions. 

The angular momentum analogue of the product of two forms, say g(j), can be 
found from the multiplication rule for null spinors, (19). Thus, 



(58) 



The statement is that ip, = gcf), is apolar to / if ^ is, i.e. , that 

ii J2 + is (ii ~ ^2 ~ is) 

mj^ 7714 * 

i2 is ™4 V il ^2 + is O'l - 32 - h) 

^2 "7,3 h + k )\mim^ * 



= /"^i^;™ 



is true if (57) holds. To show this, the 3j symbols arc rccoupled, using 6j symbols, 
so that j\ is coupled to J2 ^-i^d J3 to j\ — J2 ~ k- Without writing out the full 
expression, these couplings quickly show that the smallest intermediate angular 
momentum is j\ — J2 cind the largest one is also j\ — J2- This allows (57) to come 
into play and the result follows. 

Two apolar forms of the same order (class), n, are sometimes called conjugate.'^^ 
From (57), this means that the joint invariant, 

/^^^m, = (a./3)^ (59) 

vanishes. 



'^^ The terminology seems to originate with Rosanes, [61]. Reye uses the term apolarity as do 

Grace and Young, [29]. The classical expression for the joint invariant (or lineo linear invariant) 
can be seen in many places e.g. Elliott, [39], §49, Grace and Young, [29] §177, Glenn, [31], equn 
(71), Coolidge, [60] p. 369, Rosanes, [61]. The relation between (36) and (38) is on p. 58 of Elliott. 
Following Sylvester, Elliott considers polars as examples of emanants. 
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A form, 0, of order n, is conjugate to the nth power (c/ (10)) of one of its 
factors, essentially by definition of the roots of = 0. These powers are null 
forms, on my previous definition, and their coefficients are the components of what 
I termed the principal null 2j + 1-spinors of [21]. In general there are n = 2j of 
them and any linear combination is obviously conjugate to (p. Conversely, any form 
conjugate to (f) can be written as a linear combination of the principal null forms of 
(j), essentially by a coimting argument (Rosanes, [61]). 

Since odd order forms (half-odd-integer spin) are self-conjugate, they can al- 
ways be expanded as sums of nth powers. 

An easy theorem, which has sometimes been used as the definition of apolarity, 
is that / is apolar to g if, and only if, it is conjugate to the product gcf) for all factors 
(f). Thus, conjugacy reads, using (58) for the product of forms. 



If true for all then can be removed and, setting = j2 + j^, there results 



This proves necessity. Sufficiency has been shown just above. 

A specific example, that combines these lemmas, follows from the groupings of 
the equation. 



which say, in form language, that, if a" and are conjugate n-ics, then 

(ac) a^~^ and h'^~^ are conjugate (n — l)-ics. The particular 3j symbol here is the 
same as the matrix-spinor Uj^{j) {A = ±1/2) employed by Dirac in his higher spin 
equations (see Corson, [62]). 

As another illustration, in terms of 3j symbols the necessary and sufficient 
condition for there to be at least e equal roots of the binary quantic equation 

— 0, is that the apolarity. 






(60) 




(61) 



should have solutions for ^ corresponding to the equal root. This implies 



(62) 
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In particular, as a trivial check, when e = 2j, (62) implies oc ^(^j^, a null spinor. 
At the other extreme, if e = 1 {i.e. general position), (61) is just the polynomial, 
= 0, each of the 2j roots qualifying as an 'equal' root. 

In general, (61) is a set of e simultaneous {2j + 1 — e)-ics with at least one 
common root, the equal root. If there are two common roots, there are two sets of 
e equal roots of the original quantic, and so on. 

As intimated earlier, for calculational purposes there is no need to evaluate the 
3j-symbols in (61) as they are coded in the polar, the derivative form of which, 
(53) or (55), shows that the e polynomials are just derivatives applied to the form 

Explicitly (c/ Sturm, [59]), putting in a bookkeeping null spinor, a(^j_j,-^, 

(63) 

I also remark that the purely 2-spinor standpoint, (Penrose, [63]), is equally 
as effective. 

To summarise the correspondances so far between angular momentum and bi- 
nary invariant theory: Forms of order (or class) n correspond to spin-(n/2) spinors. 
The polar of one form with respect to another corresponds to the addition of the 
angular momenta to give the minimum total value. The product of two forms does 
the same thing but to the maximum value, while the transvectant interpolates and 
gives all allowed values. 

For the ternary case things are not so complete or straightforward but they are 
geometrically more interesting. I will return to this topic later and now take up the 
theory of harmonic polynomials, for which it is pertinant. 



The connection of invariant theory and angular momentum theory is, of course, well known 
but does not seem to be specifically used, to any degree. Any necessary calculations tend to be 
done again. However, for interesting modern developments, with an algebraic-geometric slant, 
see Abdesselam and Chipalkatti, [64] , and references given there. 
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5. Harmonic projection. Poles. 



Historically, the Laplace coefficient P^(cos7) arose in the Taylor expansion 
of the 'displaced' 1/r potential, or static Green function. 



1 _ - (r/)-Pjcos7) _ H^jr^r') 

' ' n=0 n=0 

= e-^'-^ - = y (r'.V)- - 

r ^ n! ^ ^ r 

n=0 



(64) 



yielding 



if^(r,r') = ^^r2^+^(r'.V)^-, (65) 
L! r 

given by Thomson and Tait, [8] p. 157 equn.(54). Symmetry produces the equivalent 
expression, 

H^{r,v') = tn^r"^+\r.Wr^- (66) 

A form which displays the symmetry can be found by remarking that, as we 
know (7), the left-hand side of (7), {a^x + ayi/ + Ci-,z)^, is a spherical harmonic, of 
degree L. This can be checked by direct differentiation, assuming the null condition, 
a . a = 0. Hence, in particular, 

(r.V)^ i (67) 

is a spherical harmonic in r, in agreement with (66) of course. The next step involves 
an application to (67) of the useful differentiation theorem, due to Niven, [65], see 
Hobson, [10], and [12] (6.169), 

Ynir) = (-1)" r^''^' Yni^) ' > (68) 
where is any solid spherical harmonic. There results Niven's interesting form, 

I note the similarity of the equivalent expressions, (69) and (65), with the polar 
forms in (54). 



I do not bother with conditions or regions of convergence. 
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The expression of spherical harmonics as derivatives, normal or fractional, of 
1/r was taken as fundamental by Thomson and Tait, [8], and also was adopted and 
extended by Maxwell, [9]. 

In this connection, it is best to begin historically with what might be called 
Gauss' harmonic expansion. Gauss proved constructively, in a brief 'Nachlass', that 
a rational, integral polynomial function, of degree n, /„(!■), can be represented by 
a finite sum of solid spherical harmonics, F^(r), as, 

[n/2] ^7Q^ 



s=0 

By repeatedly hitting the left-hand side with the Laplacian until zero is ob- 
tained and then solving the resulting equations from the bottom up, all the har- 
monics, can be found in terms of /^.^^ The expansion of a given is, therefore, 
unique. 

To find an explicit expression for the in terms of which would be de- 
sirable, one can start by obtaining the harmonic, Y^, from This can be found 
in the work by Clebsch, [67]. The relevant theorem, in the way he states it, can be 
motivated, in our terms, by writing (70] as. 



= + +r'Y,_, + ...) ^^^^ 

'n-2 ) 



n ' ' J n 



where the first line is the iteration of the second. This formula expresses the stan- 
dard decomposition of (the space of) homogeneous polynomials and is usually de- 
rived first {e.g. Vilenkin, [33], Berger et al, [68]) with the expansion (70) as a 
consequence). 

One refers to Y^ as the harmonic projection of Y^ = H{f^), and Gauss' 
expansion can be construed as a (finite) series of harmonic projections, 

/„ = HifJ + r' H{f^_,) + r^ H{f^_,) + ... (72) 

Clebsch, [67] eqn.3, gives the harmonic projection, although he does not use 



A specific example is worked out by MacMillan, [11], §207. See also Heine, [66] vol.1, pp. 324- 
325. 



23 



this terminology, in the form, 
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J.2 ^ ^4 ^2 



rT + 



fm ■ (73) 



2(2m-l) 2.4(2m- l)(2m-3) 

Exhibiting all the projections, H{ f^), in terms of the starting out function, 
will complete the evaluation. To do this, I follow the Gauss procedure, by writing 
out the series again, 

/n = -^n + + • • • + r"^* ^n-2s + r^*"*"^ ^n-2s-2 + • • • ) 

and then acting on it with A^, which kills the first s terms. This leaves 

A = ^ H ^n-2. + 5 ^'n.-2.-2 + • • • 

where Ag{n) is the constant, 

A^(n) = 2.4. . . .2s (2n - 2s + l)(2n - 2s - 1) ... (2n - 4s + 3) (74) 

and i? is a constant that I do need because the harmonic (first) part of the equation 
gives, 

yn-2s = ^(/n-2 J = H{A^ fj . (75) 

In this way, all the terms in (72) can be considered calculated. Equation (75) 
with (74), was obtained by Prasad, [71], and later by Dougall, [72] 

Prasad goes on to express the expansion in a compact way using the more 
general differentiation theorem of Hobson, [73], [10], §80, 

on „I 1 
H{fn) = (-ir(^^'"^'/n(V)-, (76) 

which is an extension of Clebsch's result, (73). (The more specific result (68) follows 
from (76).) Easy substitution yields the final expansion, 

/„W = (-1)" E K/n(P)];. 



A more sophisticated treatment is given in Vilenkin, [33] Chap. XI. Also the quantities and 
A are reciprocal, in the sense of Salmon, and the projection operator takes a typical form. 
Compare with the Hilbert projectors, (48) (49) (52). is a raising, and A a lowering, operator, 
cf Elliott, [69]. The Lie algebra aspect of spherical harmonics is outlined by Howe, [70], as an 
example of a bigger scheme. His other examples are also instructive. 

There is an ellipsis of an ellipsis in the expression for in Hobson [10], §96. 
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where p is the vector of gradients, p = {a, P, 7) = {d^, dy, d^), and 

92 92 92 



P da^ a/32 ■ 

As noted by Clebsch, specific choices of /„, produce known harmonics. For 
example 

/Jr) = C (r.r'r , ^ = 2-2-M! ^ ll^lllll^IiJZl) (77) 

nlnl nl 

yields the usual Laplace coefficients, P^(cos7), as follows immediately from (76) and 
(65). Hobson, [10], uses (76) to dcvclopc the properties of the zonal and tesseral 
spherical harmonics following Maxwell, and says everything that needs saying. 

Now (r.r')"^ is a rather special nth order polynomial (in r) which suggests a 
natural extension to the general product of n factors, 2^ 



/n(r) = C' n (P« • ^) ' PW " (^(i)' ^(i)) (''8) 
whose harmonic projection, with C as in (77), gives the solid spherical harmonic 

>;c-.p(.,) = ^'-'"^'n(p(.)-v) 

n 



(79) 

f_-i\n n ^ \ / 

n! . ^ . 

1=1 

where Vr, is a directed derivative. 

Maxwell refers to the points where the rays, determined by the vectors P(j), 
intersect the unit 2-sphere as the poles of the harmonic, Y^, and considers the 
harmonic to depend on them. I have indicated this by the dependence on the unit 
vectors P(^iy As noted by Maxwell, Gauss, in the same Nachlass referred to earlier, 
sketchily outlines a geometrical meaning of 'Sphere functions' referring to the poles 
as 'bestimmte Punkte'. 

Conversely, by a non-rigorous counting method. Maxwell shows that any ratio- 
nal integral harmonic, of order n can be represented as a multipole, (79). This 
is an important point and is the content of Sylvester's theorem. 



2^ Usually the vectors p^j) are taken to be unit ones. 
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6. Clebsch— Sylvester theorem. 

Proceeding a little more carefully, cf [10], [73], [74], [75] §51, the question is, 
given a specific harmonic, Y^, can one find a general homogeneous polynomial, 
of product form 

n 

/n(r)=Cn(P(i)-r), (80) 

of which is the harmonic projection? 

The relation between Y^ and is the basic Gauss expansion, (71). We see 
that this does not determine uniquely because of the the final term, which can 
be chosen arbitrarily, given just Y^. It is this 'gauge' freedom that can be used to 
work into the product form, (80). Saying it again, we are looking to write the 
harmonic Y^ is the form, 

n 

yn{r)=Cl[{p^,yr)+r'G^_,{v), (81) 

i=l 

by appropriate choice of the general homogeneous polynomial of order (n— 2), G„_2- 
This is a purely algebraic problem, first raised, and solved, by Clebsch and, later 
clearly, but more cursorily, by Sylvester. The question is independent of the notion 
of poles. 

Interpreting (x, y, z) as projective (homogeneous) coordinates in the plane, the 
terms in (81), correspond to lines h^, P(j-) . r = 0, a curve, C^, of order n, Y^{r) = 0, 
a curve, C„_2, of order (n — 2), G^_2 — 0, and an absolute null conic, C^, = 
r . r = x"^ + + z"^ = 0. Geometrically (81) says that the n lines /ij meet the curve 
C„ in the 2n points where the conic C2 cuts that curve. These points occur in n 
complex conjugate pairs and it is only lines joining such pairs that are real lines, 
and so can be identified with /i^. There is therefore only one way of doing this, so 
that the decomposition into a real product is unique. This is commonly referred 
to as Sylvester's theorem but there is a very clear and explicit statement of it in 
Clebsch, [67] p. 350, and so I suggest it be called the Clebsch-Sylvester theorem. 
The problem may be that Clebsch uses a rather complicated elimination procedure. 

A further piece of old fashioned geometric terminology is to say that the line, 
h, is the polar, with respect to the null conic, of the point (the pole) whose homo- 
geneous coordinates are the line coordinates (direction cosines) of h?^ 

There seems to be a curious congruence of terms as the word pole is used both in its projective 
geometric sense and in its physical sense, although both derive from the same root meaning 
'axis'. 
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The determination of the intersections, i.e. of the poles, by ehmination, is 
a separate algebraic matter best approached by first using the Cartan map to 
parametrise, this time, the variables, {x,y,z), rather than the coefficients, of the 
ternary n-ic, Y^, 

^n(^>y>^) = ^^j^Sgra^V^""; P + q + r = n, (82) 

so that the null conic, C2, is automatically satisfied. 

In my previous notation I denoted the null r by b, with the Cartan map, or 
uniformisation in algebraic geometry. 



b™ — lb,, „ b™ + zb,, „ 

^'^~^V^^^ ' ^-'^^V^^"^ ' bo = -ib, = v^^ry, (83) 

which differs, only by conventions, from existing choices, [74], [75]. 

Before pursuing the required elimination, I wish to make some comments re- 
lated to (83). Substitution of (83) into the ternary n-ic, (82), yields a binary 2n-ic, 

^n(t'.,^,y =>^2n(e,r/) = J]^a,,^^r/^; p + q = 2n. (84) 

Equation (12) can be considered as a formal expression of this equivalence 
through the equality of the symbolic forms of the two types of quantics. In fact it 
is helpful, in the present situation, to draw attention to the similarity between the 
true power (r' . r)" of (77), whose harmonic projection gives and the symbolic 
form (a . r)"' of a general ternary quantic.^° Then, still symbolically, the harmonic 
projection yields 



^' ^ (85) 



( „2n+l vyn ^ 



nl r 



The equivalence between the invariant theory 

quantic, which is what we really have here, has 
mention of the work of Cartan on spinors nor 
Salmon. 

See Problem 9, [10], p. 176. cf [67]. 



of a binary quantic and an orthogonal ternary 

been utilised by Littlewood [76]. He makes no 
of the relevant considerations of Burnside and 
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7. Computation of the poles. 



The elimination, to find the factors, i.e. the poles, is reduced to solving the 
2n-ic, 

>^2n(e,^) = 0, (86) 

for the ratio t = ^/rj. From reality, for every root, t, there is another equal to —1/t*. 
This means that the solutions for the corresponding vector b (projectively a point) 
are complex conjugates, as mentioned above. The line, h, joining two such points 
is. 



X 


y 


z 


K 


K 


K 


K 


K 


K 







(87) 



or 



h: r.(bxb*)=0, 
so that the real pole vectors, p^^^, in (81) are, 

P(i) = ±i b X b* , 

for every pair of null solutions b, b*. The direction cosines (line coordinates) are 

b X b* 



P(i) = ±^ 



b.b' 



(88) 



The spherical harmonics, (79) now look like, 

b 



b 



idx ^dy 

^cc \ i>z 

b* b* b* 



(89) 



As an example take the real harmonic polynomial, 

H{v) = -3a;^ - - 8z^ - Gx'^y^ + 24y'^z'^ + 24x^z^ - 60V2x'^yz + 20V2y^^ . 

Replace r by b and then use the Cartan map to get the octavic (without an 8th 
power), 

H{C, v) = 20^r;(2zV2(e^ + ry^) - rv^f) = , 
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whose roots are easily obtained. I give only one, 

^/ I fs i 



2V2 2V2' 

Substitution back into the Cartan map yields the null vector, 
^_ v^^Sv^z 5^2 V6i 1 VSi 



8 8 ' 8 8 ' 2 2 
and evaluation of the direction cosines (88) gives. 



P(i) 



2 ^2 1 
3' 3~' ~3 



In fact, in this case, the poles are the corners of a regular tetrahedron, centred 
on the origin with one corner on the 2;-axis, one lying in the a;y-plane and one edge 
parallel to the x-axis.^-*^ 

The computation of the null vectors can be bypassed by going back essentially 
to the factorisations and the basic equality (84) which is. 



1=1 1=1 ^ 



(90) 



On rewriting both sides, the left one best in spherical coordinates, = = rf't^-, 
b_i = rf, bg = V^^r] = \/2rf't, one finds, 

n n 

n K Pw) = n P(<) + P(<) + PS) ) 



i=i \ II 11^ 



From this identity the values of the spherical (contrastandard) components of 
the pole vectors can be read off. When normalised, I find 

5(.. = -TTw(''.;^(i-l«,t), (92) 

where I have labelled the pole by the root, t (dropping the index 'i'). The transfor- 
mation t ^ —1/t* just reverses the direction of the ray, giving nothing new. 
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The converse constructive calculation of -^(r) is given in MacMillan, [11], p. 397. 
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For convenience, I give the Cartesian components as well, 

p(*) = r^(^^^^' -i"^^' -^(i-N'))- (93) 

The use of the Cartan map to solve the elimination problem, leading to the 
result (93), is probably the most efficient technique of finding the poles. It occurs in 
Backus, [77], and Baerheim, [78]. The articles of Katz and Weeks, [2] and Weeks, [3] 
contain a handy independent reworking of the basics and some relevant numerics. 

After several steps of elimination, not directly involving the Cartan map, Cleb- 
sch, [67] , also provides a means of finding the poles by first factorising a polynomial 
constructed from the harmonic, giving what I have denoted by b. The poles are 
then obtained from (87). 

More explicitly, in his notation, f2 is a homogeneous harmonic polynomial of 
order n, and Clebsch proves that the (finite) expression, of order 2n, 

r^il^-!i-'Ail^ + ^f±^^'^A-il^-.... (94) 
n 2! n(n-l)4! ^ ^ 

is always factorisable into 2n linear factors. These are the bs. 

As an example I choose the simple, n = 2 harmonic, Q = xy + yz + zx, and 

find, 

T =x^-2x^y-2x^z+3x'^y^ + 3x'^z'^-2xy^-2xz^ +/ - 2y^z + Sy'^z^ - 2yz^ + z^ 

= \x + y + z — xy — yz — zx) 

y + z ^ VSi\y - z\\'^ f y + z VSi\y - i 
= \ X 1 \ x 



(95) 

giving the null vector (in Cartesians), 

1 - V3i _1 + V3i 
' 2 ' 2 

and its complex conjugate, twice. Computing the poles via (88) gives -^(1, 1, 1), 
twice. For comparison, the alternative method can be seen later, in connection with 
another matter. 

In the general situation, the linear factors can be found by first setting z = 1 
and then, in crudest fashion, substituting y = mx + c into the identity, T = to 
give T(a;) = 0, Va;. In fact, one needs only two terms, say T(0) = and T'(0) = 0. 
The former determines c as the solution of an ordinary polynomial of degree 2n and 
the latter fixes m linearly in terms of c. This method is, perhaps, less efficient than 
the one based on the binary form because of the redundant information. 



There is a misprint on p. 349 in [67], n factors being stated. 
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8. Normal forms. 



In many cases, the harmonic polynomial will be presented as a combination 
of the standard set of spherical harmonics, C|^(r). This is a 'normal form' of the 
harmonic, [10] §87. Starting from this combination allows a re-expression of the 
factorisation theorem. 

For a specific order L, let the harmonic, 

$^(r) = 0f C|^(r), (96) 
be the real data under consideration. (I called this y„(r) before.) Using 

(C^(r))* = Cf(r), (97) 

reality implies the condition, 

,*M / ',\L-M i—M iL /^o^ 
(-1) = </'M , (98) 

familiar from mode decompositions into spherical harmonics in many physical situ- 
ations such as hydrodynamics. In the terminology of [21], the raising and lowering 
operator is a charge conjugation operator. 

Replacing r by a null vector, b G C"^, the real harmonic, (96), turns into, 

^M = <^¥cUh), (99) 

still with the reality condition (98) on the data. Quite generally, one has $£,(b)* = 
$^(b*), which shows that if b is a solution, i.e. if 

$^(b) = 0, 

then so is ±b*, as Clebsch and Sylvester said.^^ 

It is a useful exercise to confirm this, since the replacement r — > b, modifies 
the reality properties, (97), of the harmonics. One can use either b or (^,7/) of (83) 
to express this. We also have the explicit formulae (8) and (9) to check things out. 

Taking the complex conjugate of gives, 

$^(b)* = </>*f {cu^)r = (-i)^+^</>z^(ci,(b))* 

= (-l)^-^</>f(C^^M(b))* 

For clarity, b* here refers to the vector with complex conjugated Cartesian components. 
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is an easily found binary 2L-ic in the variables ^ and rj {cf equations (29) 
and (84) in a different notation) whose roots correspond to the Maxwell poles. 
Although the algebra and the properties of the polynomial have been given before, 
I rederive them in a not really different way. 

Looking at the polynomial equation (equivalent to (86)) as one in t = ^/ry, 
the replacement b ^ b* corresponds to charge conjugation, ^ — > —rj*, rj ^ 
and so, if t is a root, then so is — as I said before. I have thus obtained the 
factorisation displayed on the right-hand side of (90). To determine the poles, this 
2-spinor factorisation is simply rewritten in terms of vectors, as on the left-hand 
side of (90). 

The spinor ip'^ = (—?]*, ^*) is the spin conjugate spinor to tp = (^) and one can 
refer to the complex conjugate null 3-vector b* as the conjugate vector = b*.^^ 

I remark that this analysis depends on the fact that the harmonic data is real. 
I also remark that pairing other roots produces other, complex, vector factors. The 
non-triviality of the result lies in the fact that the factors are real (if $ is). Any 
binary quantic can be factorised and hence so can $(b), but with complex factors, 
in general. Some further remarks on this can be found in [82]. 

If r is reinstated, the factorisation statement, (81), is that the harmonic poly- 
nomial $i,(r) can be decomposed as 

L 

^L(r) -CHir. p(,)) + r'G^_,{r) . (100) 

i=l 

Both the polynomial, Gj^_2, and the constant, C, up to scaling of the p^j-j, are 
uniquely determined by 

The constant C can be found by substituting a value for r that lies on the conic 
and then G follows by subtraction. It can also be constructed from the poles, as 
will be shown later. 

I also make the important point that the factors are real, if is. 

For completeness, the constant C can also be determined as follows, {e.g. [75], 

It is of course also a ternary form as mentioned earlier. 

To determine the poles, it is not necessary to convert this to a polynomial in (a;, y, z) as in 
Weeks, [3]. 

"^^ A detailed discussion of the notion of spinor and null vector in connection with rotations in 
3-space can be found in Kramers' nice book on quantum mechanics, [79], based on his earlier 
papers, [80]. The definitions mean that his null vectors are i times mine. He also usefully 
extends to the relativistic case. His approach can be found more fully in Brinkman, [81]. 
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§51, [74]). Consider the pencil of curves, 

L 

^L(r)-An(r-P(.))=0, (101) 

for varying A. Each member of this pencil goes through the 2L intersections (of 
C2 and Ci). Select a point, not one of these intersections, on the conic, C2, and 
adjust A so that the corresponding curve of the pencil goes through it, which is 
possible because none of the lines r . p^^^ = vanishes at this point. Since this 
curve of the pencil now has 2L + 1 points in common with the absolute conic, C2, 
it must completely contain it, and so degenerate into the conic and a curve of the 
in — 2)nd order, C^_2- This implies that the left-hand side of (101) must have 
as a factor, the sufficiency of which is obvious. (A proof of the necessity was given 
by Ostrowski, [83], [74].) The so determined A is the sought for constant, C, the 
value being independent of the initial point selected. 

9. Apolarity and harmonic polynomials. 

This development can be given an older analytical-geometric terminology and, 
to this end, I return to invariant theory. The relevant statement is that if a ho- 
mogeneous polynomial in (a;, y, z) (a ternary form) is harmonic, it is apolar to the 
fundamental (null) conic, C2, and conversely. This is easily shown from the defini- 
tion (53) which says that the polar is obtained by replacing the line coordinates, u, 
by the derivative, V, in the envelope equation of the conic, and letting it act on the 
ternary form (curve). In my conventions, the point locus equation oi S = reads 

00001 10 O 

Sy. = X + y + z = r^r -|- r_ir~ + TqT = 2rir_i — Tq = 
and its envelope equation, S = K2, 

00001 10 O 

Up = u +v +w = + u_iU~ + UqU — 2uiU_i — Uq = 

where I have given both the Cartesian and the spherical forms. The latter is the 
preferred (canonical) form in analytical geometry^^ and what we have referred to 
as spherical coordinates are trilinear coordinates with respect to a specially chosen 

3'^ e.g. Salmon, [37], Sommerville, [55], Todd, [51] §3.3, Meyer, [84], Grace and Young, [29], Semple 
and Kneebone, [58]. 
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triangle. For my normalisations, the locus and envelope equations have exactly the 
same form.^* The replacement u — > V can be made in any coordinate system and 
we have for the polar, 

h= ^. Aag = Ar„(r) = Q, 

since is assumed harmonic. This is the announced result whose significance is 
that the construction of curves apolar to a conic is a known topic in invariant theory. 

Using suitable parametric coordinates, the analysis of such forms can be re- 
duced to that of binary forms. The equation of the curve, C^, factorises into two 
binary forms of order n, Grace and Young, [29] §§ 241-244, Schlesinger, [86]. I give 
a discussion since the material is not generally familiar. 

The parametric coordinates of a (complex) point, k, on the fundamental conic, 
S, are just the Cartan map, 

ri = K,i, rQ = V2KiK2, r_i = K2. (102) 

(I use the notation of [29] and [86]. To compare with (5), ^ = and rj = ^2-) 

If the straight line it* = cuts the conic in the points A and fx it is easily 
shown that, 

u^=X2H2, w° = (Ai/i2 + Aa/^i), w"^ = Ai/x^ . (103) 

The tangents at A and /x meet at the pole point, 

^i = Ai^i, ro = (Ai^2 + A2^i), r_i = A2//2- (104) 

Using raising and lowering, one might dispense with the u symbol, since u'^ = r'^, 
the dual with respect to the conic. 



It corresponds to the choice of fc = 2 in [58], Chap.V case 4. A treatment that uses the general 
form for the base conic is given by Lindemann, [85]. It is interesting historically to note the 
polite way these authors claim priority. This reflects the importance of invariant theory at that 
time. 

This standard way of representing a point and a line by a pole and its polar with respect to 
a fixed conic seems to date back to Hesse in the 1850's following Poncelet, and others. It was 
used by Darboux and there exists a a well known higher dimensional (spin) generalisation due 
to Clifford, involving rational norm curves. 
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Equation (104) can be written as in (16) but this time with two two-spinors, 
A and /i, 

If one requires the Cartesian coordinates to be real, then = and A2 = Hi 
which means that A = /U*, i.e. the complex conjugate transpose, which is what is 
normally referred to as the 'adjoint', or hermitian conjugate, /x^. 

Schlesingcr gives an instructive way of finding the equation of an envelope of 
class n apolar to the base conic, which, as I have said, corresponds to finding a 
harmonic homogeneous form. This condition is equivalent to the envelope being 
conjugate to all curves of order n that contain 5 as a factor. 

Let c" be such a curve, then the necessary and sufficient condition that it 
contain 5' as a part is that, from (102), 

{c^nj + c°\/2AtiK2 + c-^Avi)" = (105) 

for each Ki and The trick now is to think of the symbolic coefficients, c, as 
variables, u, and Ki and K2 as fixed coefficient quantities, in a way to be explained 
shortly. The resulting class curve (envelope) will be conjugate to by virtue of 
(105), with c ^ u, and hence will be apolar to S. 
Thus we substitute, from (103), 

= A2//2 = (V)i 
c° = --^ (A1//2 + ^2/^1) = (V)o 
= Xifii = (A/i)_i 
into (105) which factorises into two, 

(A2«;i - Ai«;2)" (a*2«^i - /^i«^2)" (106) 

with A and fx as variables. 

The important point is that, since the k are fixed, arbitrary parameters, one 
obtains an apolar curve if one replaces the different products of the k, by any quan- 
tities. In other words, Ki and K2 can he replaced by the symbolic coefficients of a 
binary form so that the above factors become (part) binary forms. Putting, pre- 
cisely, Ki = a2 and K2 = — produces the polarised binary form, 

a>- = 0, (107) 
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as the envelope equation of a curve of the nth class apolar to the fundamental conic, 
S, as soon as one sets tt* = (A//)j, using the symbolic, 



«A<= [«?(V)-i-V2aia2(A//)o + ai(A//)i]". (108) 

Precisely the same equation holds for the locus of the apolar curve, but now 
the A and fi are translated into point coordinates by the dual (104). 

When X = fj, = K, the equation, afP = 0, gives the points of contact of the 2n 
common tangents of the apolar curve and the conic. On the other hand, the locus 
equation gives the points where the curve intersects the conic. The expression in 
(107) is the nth polar with respect to A of a^". 

What this result says is that a homogeneous harmonic polynomial in x,y,z 
can be reconstituted from its evaluation on the absolute conic 5' : = 0. I will 
check this by performing a round trip from the simple harmonic polynomial, xy + 
yz + zx, through its evaluation on the absolute to give a binary quartic, thence by 
polarisation to a product of two binary quadratics and finally back to xy + yz + zx. 

For the purposes of solution, it is useful to introduce the non-homogeneous 
ratios (the Darboux coordinates), 

/Xi : = t : 1 , Ai : A2 = s : 1 

so that, e.g., 

r,:ro:r_^ = st:^{s + t):l. (109) 
On the absolute, S, s — t and (102) holds. Therefore, on S, 

xy = + - 1) ^ ^ -^(.2 ^ -^x ^ ^/^2 _ -^x 

and so 

{xy + yz + zx)\g = ^j'' + {l + i)t^ - (1 - z)t - 1 , (110) 

which is a binary quartic. Polarising this according to (107) one finds, up to a 
factor, 

t'^ s'^ - {1 - i) St {s + t) - {1 + i) {t + s) - 1 , 

and substituting (109) one regains xy + yz + zx. 

The solutions for the points on S follow by rewriting (110); 

(2C'-2C- 1)2 = 0, with ^ = (l-^)C, 
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so that 

C=^(1±V3), 

indicating two double roots which are switched under t — > —1/t*. The corresponding 
Maxwell poles are proportional to (1,1,1), as we know, which could have been 
deduced immediately from the trivial identity, 

{xy + yz + zx) + ^ (x^ + y'^ + z"^) = ^{x + y + z)"^ . 

There are three resolutions of 

{xy + yz + zx) - X{x'^ +y'^ + z"^) (111) 

into two linear factors, A being a solution of the resolving cubic, 4A^ — /A + J = , 
of the binary quartic, (110), corresponding (in this case^°) to xy + yz + zx via the 
Cartan map. / and J are the standard invariants of this quartic, here, / = 3, J = 
— 1. The real resolution occurs for A = —1/2, (twice since = 27 J^). The complex 
one, for A = 1, is exhibited in (94). 

Schlesinger calls the unique apolar curve determined by the binary form, the 
curve associated with the binary form. Thus, given an arbitrary curve of the nth 
order, — 0, there is associated with it, via its 2n intersections with the conic, 
S", a unique nth order curve, = 0, apolar to S. This gives a geometrical 

interpretation of harmonic projection and of its uniqueness. Furthermore, there 
are, in addition to the 2n points of intersection, on S, of = and = 0, 

another n{n — 2) through which it is known^-*^ that a curve, m^~^ = 0, of order 
(n — 2) can be drawn, that is, 

f^ = P<oJ;+ slmT\ (112) 

which is the equation expressing the harmonic projection f ^ aa (with (104) to 
relate x and A,//), equivalent to Gauss' formula, (71). 

A system of n lines, L", through the 2n points on S can be considered to be a 
curve of the nth order (an n-side) and could be used as a specific in (112). The 
reverse question is whether, given a form apolar to 5, one can find a such that 

alal^pLl+ slmT\ (113) 

The general quartic is discussed by Burnside, [87], who gives some geometrical interpretation. 

This is Noether's intersection theorem. Roughly speaking, a curve, F, order n, that goes 
through the intersections of curve, C, order n, and conic, S, order 2, is F = pC + MS, where 
M is a curve of order n — 2 which goes through the remaining intersections of F and C. 
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but this would seem to be obvious because the apolar ternary form, a" a J], deter- 
mines the binary form a^^ (and vice versa) which fixes the 2n points on S from 
which the n-side, L^, can be constructed. Then (112) can be employed giving 
(113).^^ Of course, this is exactly Sylvester's approach to Maxwell's poles, which is 
equivalent to Clebsch's earlier theory. 

There are many theorems regarding this set up. For example, it is easy to show 
that if two binary forms are conjugate, then the associated curve, i.e. harmonic 
ternary form, of one of them is conjugate to any curve, i.e. ternary form, passing 
through the 2n points (on S) corresponding to the other. Conversely, if, through 
2n points on S, one can draw a curve, 6^ = 0, that is conjugate to the apolar curve, 
— 0, then every nth order curve that can be drawn through these points has 
the same property and the binary form of these 2n points is conjugate to afJ^. 

For, if = is a second curve which goes through the 2n points, then 

However, both curves on the right-hand side are conjugate to = 0, the first by 

assumption and the second because aj] = is apolar to = and so, therefore, 
is c^. Thus a system of n lines containing the 2n points represents a conjugate n- 
side of a'^a'^ = which cuts the conic in 2n points whose binary form is conjugate 
to a^^. 

In terms of s and t, the conjugate transformation A^^ — > — //2) ^^2 ~^ A*i is 
t — > —1/s*. The point coordinates are real (in Cartesians) if ts* = —1 and there 
are 2 real degrees of freedom, as required for the real plane. 

For applications, one has to impose the condition that the curve be a real one, 
i.e. that the form, (108), that gives it is real. 

10. More angular momentum theory. 

After this brief foray into the realms of invariant theory and simple analytical 
geometry, I return to angular momentum theory by remarking that in this, one 
concentrates on the coefficients, 0"^, of the forms, see (33), as exemplified in my 
earlier discussion of apolarity using 3j-symbols. Writing the harmonic as (96) allows 
one, sometimes, to exploit the considerable amount of existing angular momentum 
algebra. 

This way of replacing a binary form by two ternary ones is outlined in Salmon, [37], §190. See 
also Burnside [87]. 
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In dealing with higher-spin quantities, which is what $ is, there is always the 
option of using a (multi) two-spinor description. Since two-spinors are the defining 
representation there are few basic algebraic operations, which is a big advantage, 
but calculations can become longwinded. Penrose's, [88] , treatment of the algebraic 
structure of the Riemann curvature is, however, a good example of the elegance that 
can be achieved (Pirani, [89]). A well known alternative, for integer spins, is to use 
'bivector', or 3-spinor, space, in which higher spins are represented by symmetric, 
traceless tensors, e.g. 4>m^...m-^ ("^i = I'O? —1).^'^ 

In this connection, I firstly note that the null 2j -|-l-spinor, (9), which comprises 
the variables of the binary form, can be written using 3j-symbols as, [21], 



Sim 



1 



u 



1 



(114) 



where is the rectangular matrix-spinor exhibited, up to a factor, as a 3j-symbol 
in (60) and = (^). The combination of the us is the generalised Zj symbol that 
uniquely relates the symmetrised \ ® . . .®\ representation to the spin-j one {e.g. 
Fierz, [91], Williams, [15], Ansari, [92], [12], pp.421+.]). Assuming j integral, the 
two-spinors can be combined in pairs to give the same spin-1 null vector, a, ((5), 
(16)), and, after some 3j recombinations, one encounters the familiar polynomial 
expression for the spherical harmonics. 



C^{a) = N 



where 



L nil ^ \ f ^ ~ ^ '^2 



MlL-l/V m 1 L-2 



1 nij^ 
m" 1 ol^'^i^'^^ 



(115) 



iv=n 



k=l 



4fc2 - 1 



1/2 



or, in matrix form. 



C\a) = N A^^ (L) (L - 1) . . . (1) a^^ a^^ . . . a, 
with the definition, 

L ^^^JMm [m 1 L-1 



(116) 



^"^ This is employed by Debever, Synge and others in the classification of the Riemann tensor, 
j = 2. (See [90] for comments and a {2j + 1) treatment with the binary quartic concomitant 
basis written in 3j form.) 
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This expansion is valid for any 3-vector (a 3-spinor) as it follows from the 
general recursion, 



"4L2-1" 


1/2^ 


L 





Ci(r) = 

which yields, 

C\v) = NA-^^{L) A-^^{L - 1) . . .A--(l) v^v^^ ...r^^. (118) 
However the 'inverse' formula, 

^m, • • • V, = C^a) . A^^ (L) A^^ (L - 1) . . . (1) , (119) 

is not valid for any 3-vector, r, differing by terms having as a factor. 
The complete harmonic, (96), is 

$^(r) = </>^.C^(r) 

= AT 0^ (L) A^^{L-1)... (1) r^^r^^ . . . r^^ (120) 

= (f)"^i-^L r r r 

which defines the symmetric, traceless components, of the harmonic in 

textbook fashion. It is easily checked, using angular momentum theory, that the 
product of matrix-spinors, . . . A"^^ ^ ig traceless on the ms. 

In this approach, Maxwell's poles are introduced by noting that any (2L + 1)- 
spinor can be expressed in terms of a symmetrised product of, generally different, 
3-spinor s, '(/'*, 

0^ = (L) (L - 1) . . . (1) ^l^y^^ . . . i;^^^ , (121) 
corresponding to the irreducible bivector (tensor) equation, 

(f>m,...mj = i^lm.i^m^ ■ ■ -^m^) " tracc terms . (122) 

(There are no 'trace terms' for 2-spinors.) 

The reality condition (98) can be achieved, courtesy of 3j properties, by making 
all the 'factors', '0^, obey the same requirement, i.e. they correspond to real 3- 
vectors. 
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The factors, ip'^, are the Maxwell poles, p^^^, as can be seen by reconstructing 
the harmonic (96), using (120) and (122), 

L 

$^(r) = 0^ . C^(r) = C • ^1 + ^^-2(0 , (123) 

where G arises from the trace terms and can be written out explicitly in terms of 
the -0* = p(i). 

We thus regain the decomposition structure, (100), in a slightly heavy-handed, 
but interesting, way. 

A numerical example is always helpful and I return to the harmonic polynomial 
used at the end of §5. The four poles can be found in MacMillan, [11]. They are, 
in Cartesians, 

P(i) = (0,0,1), P(2) = ^(0,2A-1) 



3 

P(3) = -^(V^,Al), P(4) = ^(V^,-V2,-l). 



Then 



135 JJ(r.p(i)) = -5z{l2V2x'^y -6x'^z-4:V2y^ -6y'^z + z^) , 

where the constant, 135, follows by substituting in values satisfying x'^ + y'^ + z'^ = 0. 
As mentioned before, the last term, G, in (123) could be found by subtraction but 
now I wish to obtain it directly from the poles via the trace terms in (122) which I 
contract with rs to get 

4 

r^i ...r^40^^...^^ = JJ(r.pW) -trace terms. (124) 

Using the standard expression for the traceless part of a fourth rank tensor*^ {e.g. 
Jaric, [95]), or from a short direct calculation, it follows that. 



n + 4: ^ (n + 4)(n + 2) 



trace terms = — ;— ;- Eg — — Eg , n = 3 , 



The tensor approach to Maxwell's poles has been discussed by Zou and Zheng, [93]. See also 
Applequist, [94]. 

'^^ Applied either to (f> or to the product of rs. A general expression is available, but I will not 
write it out. 
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where Eg and Eg are the sums of the 6 pair traces and of the 3 double pair traces, 
respectively, 

^6 = (P(3) • P(4)) + • • • , S3 = (p(i) . P(2)) (P(3) . P(4)) + • • • • 

For my example, arithmetic gives Eg = 2r^/9, E3 = 1/3 and, simply. 




agreeing with the result of subtraction. 

The fact that the remainder term in (123) can be expressed in terms of the 
poles also follows by noting that there are L{L — l)/2 constants in G]^_2 which 
could be found by substituting the L{L — l)/2 cross products, pj-j) x p^^-,, for r in 
(123), when the first term on the right-hand side vanishes. 

The standard spherical harmonic, C^j^(r), can be expressed in terms of the null 
{2j + l)-spinors constructed from the A and fj, introduced earlier, (104), 

^^SW = (-lf'l<^i±|P(:;77)AS>«, 2, = L, (125) 
with A(^) = (Ji) and = (Ml). Also, 

{a.rf^ = (X^'\^^^){-p^^H^'^)={{^2^i-^i^2) (/^2ei-/^i^2))'', (126) 
(c/ (106)). 

A two-spinor approach to Maxwell's poles and Sylvester's theorem is given by 
Dennis, [4], who makes some useful comments on the relation to other techniques. 

11. Discussion. 

Maxwell's motivation for introducing poles was to provide a coordinate inde- 
pendent geometrical characterisation of spherical harmonics in 3-space; in his own 
words (in the first edition of his book ^^) "emancipating our ideas from the thrall- 
dom of systems of coordinates". This aspect, no doubt, accounts for the recent 
rediscovery of the procedure when analysing astronomical data. 

The geometrical viewpoint can also be used advantageously when finding the 
harmonics invariant under some finite subgroup of the rotation group. Maxwell's 

It is interesting to compare the treatments of this topic in the first two editions. 
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theory of poles has been employed in this regard by Poole, [96], and Laporte, [97], 
although it is not required of course. Polyhedral harmonics had earlier been intro- 
duced by Pockels in his classic book, [98]. 

It is not clear to me whether the equivalence of apolarity and harmonicity dis- 
cussed in §9 meets Sylvester's requirement of a geometrical or algebraical definition 
of a (spherical) harmonic. In fact it is not obvious that one needs such a definition in 
view of Maxwell's rebuttal on pp. 201-202 of his second edition (without mentioning 
names) of Sylvester's objections. 

The limitation to integral polynomials can be removed. Gauss' harmonic ex- 
pansion then becomes an infinite series, as noted by Hobson, [10], p. 128. In invariant 
theory this can be covered by using Sylvester's perpetuant construction, which allows 
the order to become infinite. 
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